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Abstract 

We consider the singular boundary-value problem Am = f(u) in D, u |g_o= (f>, 
where 

1. D is a bounded C 2 -domain of R , d > 3, 

2. / : (0, oo) —7- (0, oo) is a locally Holder continuous function such that f(u) — > 
oo as u —7- at the rate u~ a , for some a £ (0, 1), 

3. and (f> is a positive continuous function satisfying certain growth assumptions. 

We show existence of solutions bounded below by a positive harmonic function, 
which are smooth in D and continuous in D. Such solutions are shown to satisfy a 
boundary Harnack principle. Probabilistic techniques are used in proving the main 
results. 
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1 Introduction 



The Boundary Harnack principle (BHP) is a key tool in obtaining many results in classical 
potential theory. Over the past three decades, there has been a lot of research on proving 
that positive harmonic functions satisfy the principle in very general domains [BB91]. 
The principle is stated precisely in Section 2. 

Our work began as an investigation into another natural question: Do solutions to the 
semilinear elliptic Dirichlet problem, 



-Am = f(u), in D, 
u = (/>, in dD. 



satisfy the BHP for any /,(/> and D C IR d ? One quickly observes that, in general sub- 
harmonic functions do not satisfy the BHP, but subharmonic functions bounded below 
by a positive harmonic function do satisfy the principle. Observing this naive way of 
generating comparable solutions we reformulated our question. 

Reformulated problem: Under what conditions on D,f, and </>, will (1) have solutions 
bounded below by a positive harmonic function h ? 

Let C(D) denote the Banach space of continuous functions on D with the supremum 
norm; C 2 (D) } the set of functions having all derivatives of order less than or equal to 2 
continuous in D. We only consider classical solutions u to (1), i.e., u £ C(D) D C 2 (D) 
such that u solves (1). A classical solution u of (1), is a positive solution if it is positive 
in D and non-negative on D. 

Under certain regularity conditions on D C IR d , d > 3 and </>, Chen, Williams and Zhao 
[CWZ94] showed the existence of solutions to (1) bounded below by a positive harmonic 
function, if / satisfies —u < f(u) < u. 

The equation Au = u' p in D, u = (j> in dD has also been widely studied. For 1 < p < 2, 
it has been studied probabilistically using the exit measure of super-Brownian motion, (a 
measure valued branching process [Dyn94]), by Dynkin, Le Gall, Kuznetsov, and others 
[[LG95, DK98]]. 

Using analytic techniques, properties of solutions when f(u) = u p } p > 1, with both 
finite and singular boundary conditions have been studied by a number of authors. We 
briefly review a sample. Bandle and Marcus [BM92] give results on asymptotic behavior 
and uniqueness of the "blow-up solution" u which includes the case of the non-linear 
problem Au = u' p for any p > 1, (Loewner and Nirenberg had studied the special case of 



p = (d + 2)/(o? — 2) [LN74]). Fabbris and Veron have studied the problem of removable 
singularities [FV96] and related work on boundary singularities can be found in [VG91]. 

In the above examples, singularities were in u and not in Am. Choi, Lazer, and Mckenna 
have studied a variety of singular boundary value problems, ([LM91], [YM98]), of the type 
Au = —p(x)u~ a in D } u = in dD } where a > and p is a non-negative function. From 
their work, existence of solutions bounded below by a positive harmonic function can be 
established. 

Attention then turned towards our singular boundary value problem Au = u~ a , u = </>. We 
cannot impose zero boundary conditions and expect positive solutions (it would contradict 
the maximum principle). What happens when we let </> be zero on part of the boundary 
and positive on the complement ? There is a wealth of literature on semi-linear partial 
differential equations. To the best of our knowledge and discussions with researchers in 
the area, the above singular problem has not been considered. 
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1.1 Problem ingredients and main result: 

We hx a harmonic function h and then proceed to list the assumptions on D } /, and </> 
under which we can produce solutions u to (1) bounded below by h . The condition on </> 
is stated as an hypothesis in Theorem 1.1. 

h - Fix a point x in D. Let h G C(D; [0, 1]) be a positive harmonic function in D } 
such that h (x ) = 1/2 and A = {x G 3D : h (x) = 0} is a connected non-empty 
subset of the boundary. 

D - We will consider a simply-connected bounded C 2 -domain D in IR d , d > 3. 

/ - Let a > and < a < 1. Let / : (0, oo) — >■ (0, oo), satisfy the following conditions: 

(Fl) if < z < a, then f(z) < z~ a ; 

(F2) / is locally Holder continuous with exponent a; 

(F3) there exists M , such that f(z) < M for z > a. 

For first reading, the reader may assume that f(u) = u~ a . 

We shall hx a, a, D } f and h for the rest of this article. Our main existence theorem is 
stated below. 



Theorem 1.1. There exists < C\ = C\(D, f,h ) < oo such that if cj) £ C(dD) and 
4>(x) > (1 + Ci)ho(x) then there exists a solution u to (1) such that u > h . 



1.2 Outline of proof 

Notation: For any x = (x 1 , x 2 , . . . , x d ) in IR d , we write x = (x, x d ), where 

x = (x 1 , . . . ,x d_1 ). Cq(D) denotes the subspace of C(D) consisting of functions which 

vanish on the boundary. 



Let (0, T) be the canonical measurable space, on which X : X [0, oo) — >■ IR is a stochasti 



1C 



process. For x £ D } let P x denote the probability measure under which X is a Brownian 
motion starting at x and E x the corresponding expectation. Let td = inf {t > : X t $. D} 
be the exit time of the path from D. 

In new results or sections of proofs, we restart the numbering of constants from ci, c 2 , . . . . 
We note, however, that the constant C\ is fixed for the entire article and will be determined 
by Proposition 3.1. 

The basic idea of the proof is the following: First, if solutions to (1) exist, then they 
can be identified with their implicit probabilistic representation (Section II. 3, page 107, 

[Bas95]) 

u(x) = E x (f>{X TD ) - E x f ° f(u(X s ))ds 

Jo 

= E x (f)(X TD ) - / G D (x,y)f(u(y))dy, 



Id 



where Gd(x } -) is the Green function of the domain D with pole at x; second, define a 
mapping T [See (8)], from an appropriate subset fi of C(D) [See (4)]; third, show that T 
is continuous, T(O) is relatively compact (Proposition 3.2), and T(fi ) C O (Proposition 
3.1); and finally, use the Leray-Schauder fixed point theory presented in [TG83], to show 
that the above implicit probabilistic representation has a solution u £ C(D). A compact- 
ness argument shows that u £ C 2 (D). The above structure for the proof was inspired in 
part by the arguments presented in [CWZ94]. 

The following three key (formal) observations are essential to the proof of Theorem 1.1. 



(01) if u > h then essentially f(u) < h^ a } 

(02) as h is vanishing on A C dD, h (x) ~ dist(x,dD), for all x near A, 

(03) the time spent by Brownian motion in a box of height 2 r is comparable to 2 2r . 



In Section 2, we discuss the boundary Harnack principle and show that the solutions to (1) 
as in Theorem 1.1 are indeed comparable (Theorem 2.1). Using (01), Section 3 provides 
the details of the above idea for Theorem 1.1, assuming Lemma 3.1 and Proposition 3.1. 
Lemma 3.1, proved via standard analytic methods in Section 4, is needed for Proposition 
3.2. In the same section, using properties of conditioned Brownian motion, we prove 
Proposition 3.1. (02) and (03) form the basis of the proof. This proposition provides the 
groundwork for generating comparable solutions and of the constant C\ arises naturally 
here. The BHP for harmonic functions and facts about conditioned Brownian motion are 
used in proving this proposition. 

Finally, in Section 5, we conclude with some remarks and a list of open problems. 



2 Boundary Harnack principle 



Suppose u and v are two positive harmonic functions on D that vanish on a subset of 
dD. The boundary Harnack principle says that u and v tend to zero at the same rate. 
Bass and Burdzy [BB91] proved the boundary Harnack principle for positive harmonic 
functions in twisted Holder domains of order a for a £ (1/2,1]. In the same article, it 
was also shown that for each a £ (0, 1/2), there exists a twisted Holder domains of order 
a for which the boundary Harnack principle fails. 

Lemma 2.1. (Boundary Harnack principle) Suppose f : IR d_1 — > IR is a Lipschitz func- 
tion with constant A > ; | f(x) |< 1 for all x £ IR d_1 ; and let 

D = {x = (S;,x d ) £ R d :| x |< l,f(x) < x d < 2} 

D t = {x £ D :| x |< 1/2, x d < 3/2}. 

There exists c\ > which depends on X but otherwise does not depend on f such that for 
all x } y £ D\ and all positive harmonic functions g } h in D which vanish continuously on 
{z £ dD : z d = f(z)} we have 

g{x) ^ g(y) 



h(x) h(y) 



Proof: cf. [Bas95]. 



The above lemma holds (with the same c\) in domains which may be obtained from D 
by scaling. When applying the boundary Harnack principle we leave it to the reader to 
find the right choice of D and D\. In [BB91], it is also shown that the principle holds 
for L-harmonic functions for uniformly elliptic operators L in divergence form. We now 
prove below a similar result for solutions of (1) produced in Theorem 1.1. 



Theorem 2.1. Let D be a bounded C 2 -domain ofE, d , and </>i, cj) 2 be non-negative functions 
in C(dD). Let u\ and u 2 be solutions of the Dirichlet problem 



-A«, = f(ui), in D, 



Ui = fa, in dD } 



for i = 1,2, which are minorized by h . Suppose O is an open set and B is a compact 
subset of O such that u\ and u 2 vanish continuously on O D dD. There exist constants 
Ci,c 2; depending only on B } and D such that 

° 2 < ^t\ < ciE Xo {^i{X TD )) for all x £ B n D. 



E X0 ((f> 2 {X TD )) u 2 (x) 

Proof: If h\(x) = E x ((f>i(X TD )) } then hi is a harmonic function with boundary data c/V 
Let Vi = Ui — hi. Now, Aui(x) = 2f(u\(x)) > for all x £ D and v\(x) = for all 
x in dD. 

The maximum principle for subharmonic functions implies that v\(x) < for all x in D. 
Using our hypothesis about u i} we have h (x) < u\(x) < h\(x) for all x in D. Similarly 
we can show that h (x) < u 2 (x) < h 2 (x) for all x in _D, where h 2 (x) = E x ((f> 2 (X TD )). 
This implies that 



Mfl < !iiM < hM. for a ii x G £,. ( 2 ) 

/i 2 (^) w 2 (x) h (x) 

Note that both /ii and h 2} vanish continuously on O D dD. By the boundary Harnack 
principle for harmonic functions, we know that there exists c 3 such that for i = 1, 2, 



hj(x) h;(xn) „ „ 

t4-t < cst^t for a11 xeBnD. (3) 



(2), (3) and the assumption that h (x ) = 1/2, imply that there exist Ci,c 2 such that 
C2 < ^t\ < ciE X0 {(j)i{X TD )) for all x £ BnD. □ 



E X0 (<f> 2 (X TD )) u 2 (x) 
The following corollary follows easily from the above theorem 



Corollary 2.1. Suppose B } } D } Ui } and u 2 are as in Theorem 2.1. There exists ci, 
depending on B } (9, D } x 0} (f)\ and (f> 2 such that 

u-i(x) u-\(y) „ ^ 

^H^ c i^7T' x,yeBHD. 
u 2 (x) u 2 (y) 



3 Existence via a fixed point argument 

Let 

n = {ueC(D)\u>h }. (4) 

By assumptions on /, for u G , there exists / G C(D) such that 

f(u(xj) < h- a (x) V l(x) for all x G D. 

We let 

A^(x) = h^ a {x) V /(a;) for all x £ D, (5) 

and 

Ca' = {h : -D — y IR : /i is Borel measurable and |/i(x)| < |A"(x)| for all x G D}. 

The following lemma and proposition play a key role in the proof of Theorem 1.1. 

Lemma 3.1. The family of functions, {Gd(x } -)K(-) : x G D} , is uniformly integrable 
over D; i.e., 



lim 



sup / G D (x,y)K(y)dy 



0, 



AcD,m(A)->0 IxeDJA 

where m is the standard Lebesgue measure on IR d . 

Proposition 3.1. There exists C\ = C\(D, f,h ) G (0, oo) such that, for u G O 

E x J? f(u(X.))ds 
sup u < 6i. (6) 

xeD h (x) 

The constant C\ in the above proposition is the constant that is mentioned in the state- 
ment of Theorem 1.1. We prove the above lemma and proposition in the next section. 
We require the following proposition about the relatively compact subsets of Cq(D). 



Proposition 3.2. The family of functions tC = {J Gx>(-, y)g{y)dy '■ 9 £ Ck} is uniformly 
bounded and equicontinuous in Co(D), and, consequently, it is relatively compact in Co(D). 

Proof: As D is a bounded C 2 -domain, we have that for all y £ D 

lim Gd{x, y) = 0. 

x—tdD 

Using Lemma 3.1, for each g £ Ck, the family of functions {Gd(x } •)<?(•)} is uniformly 
integrable. This justifies the interchange of limit and integration; hence we obtain that 

x -► / G D (x,y)g(y)dy 

is in Co(D) for each g £ Ck- Therefore each function in K, is a member of Co(D). In 
particular, x — > J Gr){x,y)K(y)dy, is in Cq(D). For each g £ Ck-, 

Therefore, the family of functions K, are uniformly bounded and converge uniformly to 
zero as x — >■ dD. For any x,z in D and g £ /C, we have by using uniform integrability of 
the family of functions {Gd(x } -)K(-)} that 



G D {x,y)g{y)dy- / G D (z,y)g(y)dy 

D J D 

< / | Gd(x, y) — Gd(z, y) | K(y)dy — > 0, as x — > z. 
Jd 

(7) 
We can conclude that the functions in K, are equicontinuous in D. □ 

Proof of Theorem 1.1: We recall that fi = {u £ C(D) \ u > h in D}. Clearly, fi is 
a closed convex sub-space of C(D). For each u £ O , dehne Tu by 



(•I'D 

Tu(x) = E x (f>(X TD ) - E x f(u(X s ))ds for all x £ D. 

Jo 



1. Note that for each u £ O , f(u(x)) < K(x) } for all x £ _D and 

£* / f(u(X s ))ds = / G D (x,y)f(u(y))dy. 
Jo Jd 



H 



ence, 



/Ci = {E { .) f ° f(u(X s ))ds : u G } C K. 
Jo 

Using Proposition 3.2, we conclude that K,\ is relatively compact in Co(D). For 
4> G C(dD) } it is well known that E x (f)(X TD ) G C(D). Therefore, we have 

for u G , T(u) G C(D), (9) 

and 

T(0 ) is relatively compact in C(D). (10) 

2. Note that Tu(x) = <f)(x) for all x G 3D. This implies that, Tu(x) > h (x) for all 
x G <9_D. By Proposition 3.1, we have for all x G D 

Tu(x) = E x <f>(X TD ) E x j TD f(u(X s ))ds 
h (x) h (x) h (x) 

> e t mx td ) c ^ 

~ h (x) 

Therefore, if <f)(x) > (1 + Ci)h (x) for all x G dD, then T(u(x)) > h (x) for all 
x G -D. We have shown that 

T(fio) C O . (11) 

3. If u n G O is such that || u n — u ||oo — ^ 0, then f(u n (x)) — > f(u(x)) for all x £ D. 
Now for u G O , f(u(x)) < ATx) for all x £ D. We have shown in Lemma 
3.1 that E x C D K(X s )ds < oo. By assumptions on _D, we have E x (td) < oo. 
An application of the Dominated Convergence Theorem implies that T(u n (x)) — > 
T(u(x)) } for all x G D and by (10), the convergence holds in the uniform norm. 
We have shown that 

T : O — y O is continuous. (12) 



We have shown in (10), (11), and (12) that T is continuous, maps O into O and that 
T(O ) is relatively compact. Theorem 10.1 in [TG83] implies that T has a fixed point in 
O . 

Therefore, there is a function u G C(D) such that 



'0 



u {x) = E X 4>{X TD ) - E x f(u (X s ))ds. (13) 



Since (f) is continuous and all the points on the dD are regular, equation (13) implies that 
u = (f) on dD. We also have u > h in D hence u > in D. 
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To show that u £ C 2 (D) } we use a standard compactness argument. We begin by finding 
compact subsets B^ of D such that f(u )(x) < k, for all x £ B^. Let 

,/ /u f f(uo(x)), for all x e B k 
fk M X » = { k, for all x £ ^. 

The functions fk(uo) are bounded and locally Holder continuous with exponent a. Define 



v k (x) = E x {cf>(X TD )) -E x f ° f k (u (X s ))ds. (14) 

Jo 

As fk(uo) is locally Holder continuous in D and bounded, the function v k is differentiable 
and solves the Dirichlet problem ^Avk(x) = fk(u (x)) for all x £ D and Vk = 4> on <9-D- 

Since u £ , we know from (5), that fk(uo) < A"(x) A A;, where K(x) = hQ~ a (x) V l(x). 
As /i is a harmonic function on a C 2 -domain D } we can hnd < fi < 1 and an M 
, (independent of k), such that sup xeD dist(x , dD) 2 ~^ fk(u (x)) < M. By Theorem 4.9 in 
[TG83] we know that there exists Mi, (independent of £;), such that sup xe jj | Ujt(x) |< 
s up xedD (j)(x) + Mi. 

The sequence u/; is uniformly bounded and |Au/; = f(u ), in _E? m , m > k. Using 
Corollary 4.7 in [TG83] and the diagonal method, there exists a subsequence which con- 
verges in D to a function v £ C 2 (D), satisfying ^Av(x) = f(u (x)) for all x £ D and 
v(x) = 4>(x) for all x £ dD. 

As fk(u (x)) increases to f(u (x)) } by the Dominated Convergence Theorem, (13), and 
(14), we have 



v(x) = E x ((f>(X T )) - E x f(u (X s ))ds = uq(x) for all x £ D. 

Jo 

The above argument implies that there exists u £ C(D)DC 2 (D) that solves the Dirichlet 
problem (1). □ 



4 Feed back analysis 



In this section we prove Lemma 3.1 and Proposition 3.1. In order to notationally simplify 
the proof, we assume without loss of generality the following restrictions on the set A and 
geometry of the domain D: 

9 



D 




Figure 1: Cross section of D 

(SI) We recall that A = {x £ dD : h (x) = 0} is a connected set. We assume that there 
exists a compact set K C D containing A, such that there exists a bounded C 2 
function F on IR d_1 such that 



KC\D C {x £ 



>d . d 



> r(x)}, 



A C {x £ R d :x d = T(x)}, 

and there is a y £ A such that A C {x £ 9_D : | x — y | < f }• 

(52) There exists r such that for any x £ <9_D, there exists a y £ D } such that -B(y, r ) C 
D and x £ dB(y } r ). 

(53) Define the vertical distance from a point x £ K D D to <9_D by 

£( z ) = z d _r(x). 

Note that since this is a bounded C 2 -domain, there exists a constant c, depending 
only on D } such that 

cS(x) < dist(x,dD) < S(x). 

Next, we define certain sub-domains of D and stopping times. For any y £ dD } let 
Q(y, a, r) = {x £ D : 8(x) < a, | x — y |< r}, 

<5i = <5(yo,i,i), 



<5(yo,i, ^)' 



to 



Q k = Q(y ,2 -fc ,l), for all k > 1, 
U k = {z G 0Q k : S(z) = 2" fc }, 

S k = {z G 9Qfc : | 5 — 2/0 | =1}; 

for any Borel set BcD, let t b = inf{t > : X t G - B} and T k = inf{t > : X t G U k } 
(see Figure 1). 

(54) Recall the parameter a, from the definition of /, in (Fl). We set a = 1. We assume 
that h (x) > | for all x G D\Qi. 

(55) Recall from (5) that for all x G D } A"(x) = h,Q a (x) V /(x). We assume that there 
exists a constant Ci such that for all x G Q\, K(x) < Cidist(x, dD)~ a . As / G C(D) } 
there exists a constant c 2 such that K(x) < c 2 for all x £ Ql C) D. Also note that 
A" is a strictly positive function in D. 



Proof of Lemma 3.1: Let e > be given. Fix y ^ x in D\Qi. Let Si > be such 
that B(y , S t ) U B(x , S t ) C D\Q t and B(y , S t ) n 5(x , ^) = 0. 

Let Ci > 0. Choose m so that Gr)(z } y ) < C\ and Gd(z } x ) < ci, for all z G Q m . 
Let z G Q m , $ > 0, and A be a measurable subset of D. We consider three regions 
Rj,j = 1,2,3, where 



i?i = B(z,8 Adist(z,dD)/2), 

R 2 = A\ (B(y , S t ) U B(z, S A dist(*, dD)/2) , 
and 

R 3 = A\ (B(x , S t ) U B(z, S A dist(*, dD)/2)) 

Let 5 = $ A dist(A, <9-D). Using assumption (S5), we have 



G D (z,y)K(y)dy < c 2 / | 2 - y | 2 d K{y)dy 

R x JB(z,s) 

< c 3 f r 2 - d (s - r)- a r d - l dr 
Jo 

< c 4 (s) 1 - Q 

< c 5 S l ~ a . (15) 
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We now proceed to estimate the integral over R 2 . Applying the boundary Harnack prin- 
ciple in D\(B(xo, S^ U B(z, 1/2[S A dist(*, dD)/2])) and in D\(B(y , Si) U B(z, 1/2[S A 
d\st(z } dD)/2])) } there exists c 6 such that, for all 
y G D\(B(x , Si) U B(z, S A dist(*, dD)/2)), 

Gd(z,v) . G D (z,y ) 

< cq-^— -, (16) 



G D (x 0} y) G D (x 0} y ) 

and for all y G D\(B(y , S t ) U B(z, S A dist(*, dD)/2j) 



G D (z } y) G D (z } x ) 

^ c e^TT— T- ( 17 ) 



G D (yo,y) G D (x ,y ) 
Hence, 



G D (z,y)K(y)dy < c 6 D ' ° / G D (x ,y)K(y)dy. 
r 2 G D {x ,y ) J R2 

Since Gr>(x 0} y) is a bounded continuous functions in _D fl -B(x ,^i) c , there exists C7 such 
that 



G D (x ,y)K(y)dy <c 7 K{y)dy. (18) 

fi 2 -/A 



Using a similar argument, we may check that 



G D (z,y)K(y)dy<c 8 K(y)dy. (19) 

R 3 J A 



As Ac RiUR 2 U R 3 , for z G Q m , (15), (18) and (19) yield 



G D (z,y)K(y)dy < / G D (z,y)K(y)dy 

A JR X 



+ / G D (z } y)K(y)dy+ / G D (z,y)K(y)dy 

J Ro J Rr 



< csS 1 -" + c 8 I K(y)dy. (20) 

J A 

For z G D\Q m} choose 
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1. 8 2 > such that B(z, 8 2 A dist(z, 3D)) n Q m -i = 0. 

2. £1,2/1 G Q m so that there exists $3 > 0, such that B(yi } 8 3 ) U B(xi } 8 3 ) C Q m . This 
will ensure that for all x G D\Q m} there exists c 9 > such that 

Gd(xi 7 x) < c 9 and Gd(vi } x) < c 9 . 

We split A into three regions, B 17 j = 1,2,3, where 

B t = B(z,S 2 Adist(z,dD)/2), 

B 2 = A\(B( yi , 8 3 ) U B(z, 8 2 A dist(z, 3D) 12)), 
and 

B 3 = A\(B( Xl , 8 3 ) U B(z, 8 2 A dist(z, 3D) 12)). 

Because K is bounded in B(z, 8 2 A dist(z, 3D)/2) } we have 



G D (x } y)K(y)dy < c w I G D (x } y)dy 

Bx Jb( x ,s 2 ) 

= c^l (21) 

Using the boundary Harnack principle again as in (16) and (17), we obtain 



G D (x,y)K(y)dy<ci2 K(y)dy. (22) 

B 2 J A 

A similar estimate, (as in (22)), holds for B 3 . Since A C B\ U B 2 U B 3} for z G Q m , we 
apply (21) and (22) to deduce that 

G D (z,y)K(y)dy < / G D (z,y)K(y)dy + / G D (z,y)K(y)dy 

A J B x JB 2 



+ / G D (z,y)K(y)dy 
Jb 3 

< c u 5 2 2 + C13 / K{y)dy. (23) 



./a 
From (20) and (23), there exists 84 > such that 

sup / G D (x,y)K(y)dy < c 14 8 x ~ a + / K(y)dy 
xED J a L 

13 



for all 8 < 84. 

For w G A, by assumption (S2), there exists Wj G D and r , such that B(w\, r ) C -D and 
itf G dB(wi 7 r ). By assumption (S5), we have 



pro 

K{y)dy < c\ I (r — r)~ a r d ~ 1 dr < 00. 

B(wi,r ) JO 

We can cover dQ\ by hnitely many balls Ni, such that each Ni touches the boundary like 
B(wi, ri) and K is bounded in ((J 8 - iV 8 ') c f] D. This implies that K G L l (D). Therefore, if 
we choose 8 small enough and A satisfying m(A) < 8 } then 



sup / Gr)(x } y)K(y)dy < e. 
xED J a 

Therefore the family of functions {Gd(x } -)K(-) '■ x G D} is uniformly 

integrable. □ 



4.1 Killed Brownian motion 



In this section, we prove some results on /i-path transforms. Consider the box Q\ in the 
domain. We denote by P® 1 the measure under which X t is a Brownian motion killed on 
exiting the box Qi and E® 1 will denote the corresponding expectation. 

If h is a positive harmonic function in Q\ then 



is the transition function of a Markov process Xh '■ X [0, 00) — > IR d U {5~}, called an 
/^-transform, or conditioned Brownian motion. We will use P^ and E^ to denote the 
corresponding probability measure and its expectation. By convention, h is taken to 
vanish at 8. 

Let Vo = inf{t > : X t G [J k £4}, W{ be the integer n such that Xy t G U n and Vi + \ = 
inf{t > Vi : X t G (J fc £4 — £V 8 }- For each Nj = Qj + \\Qj-\, we construct a subsequence 
of 14 as follows. Let m\(j) = inf{n : Xy„ G Uj}, 

n i(j) = inf{n > rrii(j) : Xy„ G ££j+i U Uj-\} for all z > 1, and 
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m i{j) = inf{n > rii-i : Xy„ G Uj} for all i > 2. 

Finally, we define a family of stopping times for each strip N 3 , namely C\ = V mt (j), 
Dj = V nt{j) . 

The following lemma concerns: (i) the time spent by conditioned Brownian motion in 
various strips of the box Q\, and (ii) the probabilities of repeated excursions to a given 
strip. 

Lemma 4.1. Let h be a positive harmonic function in Q r that vanishes continuously on 
the dQ r fl D and suppose there exists c > such that h(x) > c for all x G 3D. Then 
for x G Q r and j < r , 

1. There exist c\ < oo and p < 1 such that P^{C\ < oo) < C\p l . 

2. There exists c> such that E h (B' J - - C] I C] < oo) < c2 2j . 

tL \ i III / 

Proof: For each j, let J 3 = Q(yo, 3 • 2 - - 7 ' -1 , f ). Define 



Kj = {x G J 3 \Qj : \x-y \ >1- 2"- ? - 1 }. 

Note that Kj is the union of two compact subsets Kj and K? . Let x 8 - G dlC-DUj, for i = 
1,2. 

For each x G Uj C\ Kj, 

P x {TfJD < Tj-i) = --" > C 2 P^(T dD < Tj-i). 

As x G Uj fl K^, we can insert a box, strictly contained in Qj-\, with x at its center. The 
top and bottom of the box coincide with the Uj-\ and dD, respectively. The sides of the 
box are formed by the points {y G Qj-i :\ y — x \= 2~ 1 ~ 1 }. 

We have, P^^Tan < Tj_i) > c 3 for all x G Uj fl A"? Therefore, 



P^(r 9D < Tj-i) > c 4 for all x e Uj C\ K] . (24) 

A similar argument will show that (24) holds for x 1 a for i = 1,2. If x G A" 1 , x — > 
A^? 1 (/i(A TD ) : td < Tj_i) and x — > /i(x) are both harmonic functions that vanish contin- 
uously on dKj fl dQ\. By the boundary Harnack principle for harmonic functions, there 
exists c 5 such that 
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Px(T dD < Tj-i) > C 5 P^(T dD < Tj-i) > C 5 C 4 . (25) 

3 

Now, c 5 does not depend on j, because K\ can be obtained from K\ for n ^ /, by scaling. 
A similar estimate holds for x £ Kj. We have shown that for any x £ Uj, there exists 
c 6 > such that 



PxijdD < Tj-i) > c 6 . (26) 

Therefore, 

P*(Ci = oo) > P^(r j+1 < r,-.!, r dD < Tj ) >c 2 6 = c 7 >0 

for x £ £/j. Hence the part (1) of the lemma follows from repeated application of the 
strong Markov property. 

Given C\ < oo, the process {X t : t £ [Cj, .£)'•]} is a conditioned Brownian motion in the 
strip Qj-\\Qj + \. Following the argument on page 200, Theorem 3.2 of [Bas95], we have 
that 

E h x {D) - C) | C) < oo) < c 8 2 2j . D 

Lemma 4.2. Let h be a positive harmonic function in Q\ that vanishes on dQ r D D. For 
x £ Q r ,r < 0, 



Et 



f Ql K(X s )ds] < Cl 2 (2 " a)r . 



Proof: We have 

-OO OO / ~]j3 

V "' L C?<c 



/ r T Qi \ ~°° °° / r D i 

E h x / K{X.)ds)= £ J>£ / K(X s )ds;l C 3. 



For all t £ [C-,Dj], X t £ Q 3+ i\Q 3 _i. By assumption (S5) , it follows that K(X t ) < 
c 2 2~ ai for all t £ [C^Dj]. Hence 



E h J l K(X.)ds )<c 2 Y J 2"- ]T i#(^ - Cj; 1 £ 

^° / j=r-l . = 1 



< C 3 J] 2( 2 -^' ]T^' = C^ 2 - a >.U 



j=r—l i=l 
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Lemma 4.3. There exists c\ £ (0, oo) such that for all x\, x 2 £ Q\ H D, 

< 



EH I K(X.)ds< Cl , 

Jo 



where E x x 2 refers to the Brownian motion conditioned to stay in D , starting at X\ and 
killed at x 2 . Here ( denotes the life-time of the conditioned Brownian motion. 



Proof: Let Gd denotes the Green function of the domain D. From Theorem 5.10 in 
[CZ95] we have 



T?*2 { is(v\A f v( \ G D{x u y)G D {y,x 2 ) , 

E x\ / H^s)ds= \ h{y) -— dy. (27) 

Jo Jd G d (x 1 ,x 2 ) 

We divide the domain D into four regions: 



Rj = <i z £ D : | z — x 3 | < rj, \ z — x 3 _j \ > J> , j = 1, 2, 



where r 3 = dist(xj, dD) for j = 1, 2; 



R 3 = lzeQ 1 C\D:ze D\(R 1 U R 2 ), dist(z, 3D) < r 3 = mm C r o, r t , r 2> 



where r is from assumption (S2); and 

R 4 = D\(Ri U R 2 U R 3 ) 



We will now estimate the integral on the right-hand side of (27) over the four regions. 
Let G(x, y) = c 2 \ x — y \ 2 ~ d denote the Green function on R d , d > 3. 

By assumption (S5) , at the beginning of this section, K(x) < c 3 dist(x, 3y)~ a for all 
x £ R\. This and the 3-G inequality (Theorem 6.5 in [CZ95]), imply that 
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Ri 



G D (x 1} y)G D (y } x 2 ) 

My) — tt- ( x — d v 

G D (x 1} x 2 ) 
< cJ di8t(y,dD)-» GD{x *>? GD{ »' X2) dy 



G D (x 1 ,x 2 ) 



< c 4 \ xx - x 2 \ d\st(y } dD) a \ x x - y \ \x 2 -y\ dy 

J R x 



< c 5 / r 2 - d r d -\rx -r)- a dr 



o 



< ce s-°ds + c 7 s 1 -°ds 

Jo Jo 

< oo. (28) 



By replacing x 1} rx with x 2} r 2 in the above calculation, we obtain an identical estimate 
for R 2 . Since R 3 is a compact subset of IR d , we can cover it with balls Bi, B 2 , . . . , B n of 
radius r 3 , each of which touches the boundary of D. There exists c 8 such that Go{xi^ y) < 
c 8 , i = 1,2 for all y G U"=i Bj- A calculation analogous to that in (28) ensures that the 
integral over the region R 3 is finite. We explicate below. 



G D (x u y)G D (y,x 2 ) 

K (y) — yr- ( x — d y 

r 3 G D (x 1 ,x 2 ) 



< c 3 l di8t(y,dD)-° GD{x *'? GD{y : X2) dy 



■3 



G D (xi,x 2 ) 



< c w I xi - x 2 I / dist (y } dD)~°dy 

Jr 3 

< oo. (29) 



In the region i? 4 , K(-) } Gd(xi } ■) } Gd(x 2} ■) are all bounded continuous functions. Using 
the fact that D is bounded, it is easily seen that the integral over the region i? 4 is finite. 

□ 

Proof of Proposition 3.1 

We recall from (5) that, for all u G O , f(u) < K. Therefore, it is enough to show that 



E x ft* K{X„)ds ^ n 
sup FJ7) - Cl - 
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Let x G Q r +i\Q r . Consider a Brownian path starting at x. There are three possibilities: 

(a) the path exits the domain at its exit time from the box Q i} 

(b) the path exits the domain via dD D dQ\ but exits Q\ before that time, before, 

(c) the path exits the domain via dD D (dQ\) c . 

Therefore, we have 

E x I ° K(X s )ds = E x (f D K(X s )ds;r D = T Ql 

+ E X ([ D K(X s )ds; t d > r Ql , x TD edDn dQ l 

+E X ([ D K(X s )ds;X TD edDn (dQ 1 ) 

We analyze each case separately. 

Casel: Consider the case when {jd = tq x }. In this case, we are dealing with a Brownian 
motion conditioned to exit Q\ via dD. Hence, 



E x f / K(X s )ds; t d = t Qi ) = E x ( K(X s )ds; t d < r r+1 ,r D = r Ql 

+ ^2 Ex ( / K(X s )ds;Tj <T D < Tj +1 ,T D = t Qi 

j=r + l VO 

+ Yl Ex [ 1 (^<^Qi) Ex -j I K(X s )ds;T D < T j+1 ,T D = t Qi 

j=r+l ^ 
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If we define hj(x) = P x (tj < tqj), then the last displayed quantity is less than or equal to 



E x f / K(X s )ds;T D < r r+1 ,r D = t Qi 



(30) 



+ E 



j=r + l L J0 



E h j \ K(X s )ds 



h 3 (x) (31) 



+ ^2 E x ( 1 (Tj<TQi) E Xtj j K(X s )ds;T D < t j+u t d = t Qi ) . (32) 



j 

j=r+l ^ J0 



Note that for t < tq 1 , X t has same distribution under P x and P® 1 . As discussed in 
Section 4.1, the process X t under P® 1 is a Brownian motion killed on the boundary of 
Q\. Applying Lemma 4.2 appropriately on each summand of (31), (32) and (30), we have 
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E x ( I K(X s )ds-T D = r Ql \ < c^ 2 ~ a > 

-l -l 

+c 2 ]T 2( 2 -^'^(x) + c 3 ]T ^ (2 ~ a)3 h 3 {x). (33) 



j=r + l j=r + l 



For x G Qr+i\Qr, by definition of the levels Uj and the strong Markov property we have, 
hj( x ) — c 4% r ~ 3 ■ Therefore, we have 



- 1 



E x ( / K(X s )ds- r D = t Qi ) < Cl 2( 2 -^ + c 5 ]T 2 {1 - a)3 2 r 

j=r + l 

< c 4 2 {2 - a)r + c 5 2 r ■ (34) 



7 j=r + l 



Case2: Consider the event {tq x < T£> } X TD G 3D D dQ\}, i.e., the path exits Q\ before 
the leaving the domain D via dD D 8Q\. Then, 



E J f ° K(X s )ds; t d > r Ql , x TD edDn dQ l 
< Y E Af Ql K ( x s)ds; r D > r Ql , X TQi G sA (35) 

+ J2J2 E * ( / ^ K ( X s)ds;r D > t Qi1 X TQi G S 31 X Lq1 G sA (36) 

j k V T Qi 1 ) 

■>T D 



J k \ jL Qt 



+ J2J2 E *[ i K ( x >) ds ; t d > t Qi , x TQi g s h x Lql g s k , x TD edDndQA. 

(37) 



We analyze (35), (36), and (37) separately. 

Analysis of (35): 

Consider a summand in (35) for which j is less than r. In this j-th summand the interval 
[0, tqJ is equal to [0, Tj] U [tj, L h] _i] U [Ljj_i, tqJ, where £?j-i is the last exit time of the 
set (Qj\Qj-i) c n Qi. It follows that 
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E x lj ' K{X s )ds- t d > t Qi ;X TQi G S 3 



E x 



+ 



E m j 



K(X s )ds;T D > t Qi7 X TQi G Sj 



J 3,3-^ 



K(X s )ds + E™> / K(X s )ds 

L 3,3-1 



rrij(x 



j v^y ? 



(38) 



where m J (x) = P x (t d > t Qi ,X TQi G Sj). 
Using (34) with dD replaced by Uj, we have 

<- T 3 

E x ( / K(X s )ds;r D < t Qi ;X TQi G Sj 



< E„ 



K(X s )ds;T Ql > Tj 



< c 6 2 {2 - a)3 + c 7 2 3 . 



(39) 



The process {X t ,t G [tj, Ljj_i]}, can be viewed as Brownian motion conditioned to 
converge to a random point on Uj U Uj-\. Lemma 4.3 yields 



ET J I 3 ' 3 ^ K(X s )ds < c 8 . 

'3 



(40) 



To estimate E x 3 J ' Ql K(X s )ds, we observe the process moves inside Qj\Qj_i condi- 
tioned to exit the strip via the sides. We observe that the function is bounded above by 
a constant multiple of 2~ a ' J . Just as in Lemma 4.1, the time taken by the Brownian path 
is of the order 2 2j . We have 



E?* I Ql K(X s )ds < c 9 2 {2 -^ 3 < c w . 

L 3,3-1 



(41) 



(39), (40), and (41) imply that 



J2 Ex ( l K (. X *) ds 'i T D > T Qi , X r Ql e S 3 



'J-r 

— oo 



< 



]T c 6 2( 2 -^' + c 7 T + c 11? 



3=r 



< c 6 2^-^ r + c 7 2 r + c n ]T 



'J-r 



(42) 
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For j > r, we split the j-th term in (35) the same way as in the case j < r (38). The 
analysis is the same except for the first term. We illustrate below, the difference. If j > r, 



E x [ j ' K(X s )ds;r D >T Ql ,X TQi e SA < E a J t ' K(X s )ds) g 3 (x) 



where g 3 (x) = P x (tq 1 > Tj). Note that, g 3 (x) < Ci 2 2 r 3 . By Lemma 4.2, we conclude that 



J^eJT K(X s )ds;r D > r Ql ,X TQi G sA < c 13 ^2( 2 -^'2 r -- ? ' < c 13 2 r . (43) 

j=r ^ ' j=r 

(42) and (43) imply that 



J2 Ex ( l K (. X *) ds 'i T D > T Qi , X r Ql G Sj 



< c 6 2 {2 ~ a > + c 7 2 r + c n ]T ™j(x) + c 13 2 r + c 13 ]T m,(x) 

j =r j =r 

< c^ 2 ~ a > + c 14 2 r + c 15 P x {t d > r Ql ). (44) 

Analysis of (36): Each summand of (36) is equal to 

E x l ' K{X s )ds-T D > t Qi ,X TQi G Sj,X Lqt G S k ] 



Oi 



L Q\ 



E% / #(*,)<** «,-(*), 



Oi / 



where n 3 (x) = P x [w > tq 1 ,X t G Sj,Xl c G S^ I . Using Lemma 4.3, we have 



J Q? 



E n J \ j K(X s )ds I < c 16 . 



Hence (36) is less than or equal to 



j k 



7j2J2 Px ( TD> T ^ ' X -Qi G S i> Xl q\ G Sk ) - c ^ P ^ t d > T Qi)- ( 45 ) 
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Analysis of (37): The summands in (37) can be represented as 

E x (e Xlqi QH K(x s )ds-r D = r Q ^j -x TQi e S 31 x Lqi g s k 

From Case 1 we know that there exists a constant c\g such that, 

Ex LqC ( / K(X s )ds; t d = t Qi \ < c 18 . 

Therefore (37) is less than or equal to 






cig l^l^Px (x TQi g S^x Lqi G S k J = c 19 P x (t d > r Ql ). (46) 

j k 



By (44), (45), and (46) we conclude that 

E J f ° K(X s )ds; t d > r Ql , x TD edDn dQ l 



< c 6 2 {2 ~ a > + c n 2 r + c 20 P x (r D > r Ql ) (47) 

Case3: Consider the last case when the path leaves the domain via dD D (dQ\) c . Then 



E x ( j ° K(X s )ds; r D > r Ql ,X TD eDC\ (dQ 1 ) 



< 

3 



J2 E x ^ TQl K(X.)ds; r D > r Ql , X TQi G S^j (48) 

+ J2J2 E -\ I ^ K ( X s)ds;r D > t Qi1 X TQi G S 3 -X Lqi G Sk) (49) 

j k V T Qi 1 ) 

+ Y, E A I ° K ( x s)ds; r D > r Ql , X Lqc g S k ) (50) 



Using the same strategy as in Case 2, we can conclude (48) + (49) is less than or equal to 

c 21 [2^- a y r + 2 r + P x (r D >T Ql )}. 

Observe that in all the summands of (50), the function K(X S ) is bounded. Therefore, 
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]T E x ( I K(X s )ds : r D > TQl ; X Lq1 £ ^ J 

c 22 Yl E * \ Ex l Q c ( t d); r ^ > T Qi ; x ^ G & 



< „ 

< c 23 P x (td > t Qi ). (51) 

Therefore 

E x ( f ° K(X s )ds;r D > r Ql ,X TD £ (A)' 

< c 24 [2( 2 "") r + 2 r + P,(r D > t Qi )] . (52) 
From (34), (47), and (52) we have for x £ Q r +i\Qr, 

£* / ° K(X s )ds < c 25 [P x (r D >t Qi ) + 2^ r + 2 r ] (53) 



The function h is harmonic and vanishes on the set A. By our assumptions on D } there 
exists c 26 such that for x £ Q r+ i\Q r} c 26 2 r_1 < h (x). This and the boundary Harnack 
principle for harmonic functions imply that, for any x £ Qi, 



E x ^ D K(X s )ds ^ 
h (x) ~ 

Now for xGflfl Qj_, there exists c 28 > such that h (x) > c 28 . Since 

2 

E x I ° K(X s )ds £ C (D), 



we deduce that there exists C\ = C\(D, /, /lo) such that for all x £ _D, 

^f TD aijq)^ 



< Ci < 00. □ (54) 



h (x) 

5 Remarks 



1. The analogous results (Theorem 1.1 and Theorem 2.1) can be obtained for D C 
IR d : d = 1, 2, with minor changes in the potential theoretic arguments. 

24 



2. Boundary conditions: It is worthwhile to note that, even though the value of </> 
needs to large at certain point of the boundary there is maneuverability in the zero 
set for (f). A larger zero set for </> is accomplished by assuming a larger zero set of 
h . 

3. Open problems: We conclude this paper with some remarks on results which we 
have not obtained. 

(a) Dirichlet problem: We hx a harmonic function h and if cj) £ C(dD) and 
4>(x) > (1 + Ci)h (x) } then we guarantee existence. The more basic question 
remains unanswered: Given a </> £ C(dD) can we find a harmonic function h\, 
such that there exist solutions to (1), bounded below by hi ? Our attempts 
have so far been futile. 

(b) Critical Value of a: Is a = 1 a critical value for existence of solutions to (1) 
? We believe so, but do not have a rigorous proof. 

(c) Boundary Behavior: We have shown here that there are solutions that are 
comparable, but the broader question remains open. Do the solutions to (1) 
satisfy the boundary Harnack principle ? 

(d) Uniqueness of solutions: Are solutions to (1) unique ? Since our / is not 
monotonically increasing, there is no obvious way to show uniqueness. The 
problem seems very hard. 

(e) Lipschitz domains D: We confined ourselves to bounded C 2 domains D. A 
natural question would be: Does Theorem 1.1 hold for Lipschitz domains D ? 
Our initial calculations indicate that there is a critical Lipschitz constant. We 
shall try to make this rigorous in future work. 
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